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Abstract. In this paper, we study the requirements on orbit de- 
termination compatible with operation of next generation space 
clocks at their expected uncertainty. Using the ACES (Atomic 
Clock Ensemble in Space) mission as an example, we develop 
a relativistic model for time and frequency transfer to investi- 
gate the effects of orbit determination errors. We show that, for 
the considered orbit error models, the required uncertainty goal 
can be reached with relatively modest constraints on the or- 
bit determination of the space clock, significantly less stringent 
. \ than expected from "naive" estimates. Our results are generic 
j^ 4, to all space clocks and represent a significant step towards the 
' generalized use of next generation space clocks in fundamental 
■ physics, geodesy, and time/frequency metrology. 
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1. Introduction 

Over the last decade of the 20 lh century and the first 
few years of the 21 st , the uncertainty of atomic clocks 
has decreased by over two orders of magnitude, pass- 



ing from the 
frequency 



low 10 14 to below 
dBize S. et al. 20051 



10 16 , in relative 



Heavner T P. et al. 20051 



Oskav W.H. et al. 2006t iRosenband T. et al. 20081) . This rapid 
evolution is essentially due to recent technological break- 
throughs (laser cooling and trapping of atoms and ions), which 
allow very effective control and reduction of the motion of the 
atoms and correspondingly long interrogation times. Atomic 
fountain microwave clocks use freely falling laser cooled 
atoms and were the first to reach uncertainties below 10~ 15 
some of them being at present in the low 10~ 16 range. They 
are now closely followed, and even surpassed, by trapped 
ion and neutral atom optical clocks, the best of which show 
uncertainties below 10~ 16 . 

Space applications in fundamental physics, geodesy, time 
and frequency metrology, navigation, etc... are among the most 
promising for this new generation of clocks. Onboard terres- 
trial or solar system satellites, their exceptional frequency sta- 
bility and accuracy make them a prime tool to test the fun- 
damental laws of nature, and to study the Earth's and solar 
system gravitational potential and its evolution. In the longer 
term, they are likely to provide the primary time reference 
for the Earth, as clocks on the ground will be subject to the 
less accurate knowledge of the geopotential on the surface 
dWolf P. and Petit G. 1995h . 

For example, when comparing a clock on a low Earth or- 
biting satellite (w 1000 km altitude) to one on the ground they 
display a difference of « 10~ 10 in relative frequency due to the 
relativistic gravitational frequency shift. Measuring that differ- 
ence with 10~ 17 uncertainty would allow a test of the gravita- 
tional frequency shift to a few parts in 10 7 or equivalently, a 
determination of the potential difference between the clocks at 
the 10 cm level. The latter would contribute significantly to the 
knowledge of the geopotential and related applications in geo- 
physics, representing the first realisation of r elativistic geodesy 
i lBierhammar A. 19851: ISoffel M. etal. 19981) where the funda- 
mental observable is directly the gravitational potential via the 
relativistic redshift. 



From the above it is obvious that next generation 
space clocks at the envisaged uncertainty level require 
a fully relativistic analysis and modelling, not only of 
the clocks (in space and on the ground) but also of 
the time/fre quency transfer method used to compare 
them dAllan D.W. and Ashbv N. 19861: iKlionerS.A. 1992 : 



Blanchet L. etal. 2001 



Petit G. and Wolf P. 1994 |Wo lf P. and Petit G. 1995, 

Linet B. and T ev s s andi er P. 2002b . 



Indeed, a highly accurate space clock is of limited use unless 
it can be compared to ground clocks using a method that does 
not degrade the overall uncertainty, and unless the behaviour 
of the clocks as a function of their positions and velocities can 
be modelled with sufficient accuracy. As an example, simple 
order of magnitude estimates of the relativistic gravitational 
frequency shift show that an 1 m error on the position of the 
clocks leads to an error of ss 10~ 16 in the determination of 
their frequency difference. Similarly when using an one-way 
system (GPS like) for the time transfer an 1 m position error 
leads to an error of m 3 x 10~ 9 s in the synchronisation ie. 
> 10~ 14 in relative frequency over one day. 

In this paper, we study in more detail the requirements on 
orbit determination compatible with operation of next gener- 
ation space clocks at the required uncertainty, and based on a 
completely relativistic model. We use the example of the ACES 
(Atomic Clock Ensemble in Space) mission, an ESA-CNES 
project to be installed onboard the ISS (International Space Sta- 
tion) in 2013. It consists of two atomic clocks and a two-way 
time transfer system (microwave link, MWL) with an overall 
uncertainty goal of 1 part in 10 16 after ten day integration (see 
section[3]for more details). We show that the required accuracy 
goal can be reached with relatively modest constraints on the 
orbit determination of the space clock of w 10 m in position (for 
the considered orbit error models), which is about an order of 
magnitude less stringent than expected from "naive" estimates 
(w 1 m, see above). This is due to first order cancellation be- 
tween the velocity and position part of the orbit determination 
error in the determination of the relativistic frequency shift of 
the space clock, and to the use of a two-way time transfer sys- 
tem (MWL) which leads to first order cancellation of the posi- 
tion errors in the clock comparison (see section|6]l. Our results 
are generic to all space clocks (not limited to the ACES mis- 
sion) and represent a significant step towards the generalised 
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use of next generation space clocks in fundamental physics, 
geodesy, and time/frequency metrology, as they show that the 
constraints on the orbit determination of the space clock are 
significantly less stringent than previously thought. 

In sections |2] and [3] we briefly describe the ACES mission 
and the relativistic model used for the clocks and the time trans- 
fer. We explicitly derive the effect of orbit errors on the clock 
comparison in section [4] Up to this point, our results are com- 
pletely general (within the specified approximations) with no 
assumptions on the expected orbit determination errors. We 
then apply those results using two examples of expected or- 
bit errors (section [5]). Our main results are the calculation of 
the effect of such orbit errors on the determination of the rel- 
ativistic frequency shift of the clocks and on the time transfer 
(MWL) for the ACES mission (section[6]). We provide the over- 
all requirements on orbit determination for the ACES mission, 
and show that the mission objectives can be achieved with rela- 
tively modest performance on orbit determination. We discuss 
those results and conclude in section|7] 

2. The ACES mission 

The ACES project led by CNES and ESA aims at setting up on 
the ISS several highly stable clocks in 2013. The ACES pay- 
load includes two clocks, a hydrogen maser (SHM developed 
by TEMEX) and a cold atom clock PHARAO (developed by 
CNES) respectively for short and long term performances, and 
a microwave link for communication and time/frequency com- 
parison. The frequency stability of PHARAO onboard the ISS 
is expected to be better than 10~ 13 for one second, 3 ■ 10 -16 
over one day and 1 - 10 16 over ten days, with an accuracy goal 
of 1 ■ 10~ 16 in relative frequency. 

The ACES mission has as objectives : 

- to operate a cold atom clock in microgravity with a 100 
mHz linewidth, 

- to compare the high performances of the two atomic clocks 
in space (PHARAO and SHM) and to obtain a stability of 
3 • 10~ 16 over one day, 

- to perform time comparisons between the two space clocks 
and ground clocks, 

- to carry out tests of fundamental physics such as a gravi- 
tational redshift measurement and a test of Lorentz invari- 
ance, and to search for a possible drift of the fine structure 
constant a. 

- to perform precise measurements of the Total Electron 
Content (TEC) in the ionosphere, the tropospheric delay 
and the Newtonian potential. 

The time transfer is performed using a micro-wave two- 
way system, called Micro- Wave Link (MWL). An additional 
frequency is added in order to measure and correct the iono- 
spheric delay at the required level. It uses carriers of frequency 
13.5, 14.7 et 2.25 GHz, modulated by pseudo random codes 
respectively at 10 8 s _1 , 10 8 and 10 6 s~ l chip rates. More- 
over it has four channels that allow four ground stations to be 
compared with the ISS clock at the same time. 



According to the mission specifications, the microwave link 
has to synchronize two atomic clocks with a time stability of 
< 0.3 ps over 300 s, < 7 ps over one day, and < 23 ps over 
10 days. The performance of this link is a key issue since it 
will perform high precision time comparisons without damag- 
ing the high performances of the clocks. 

For our purposes we express the above requirements for 
the MWL in a simplified form by the temporal Allan deviation 
(c7jc(t)) expressed in seconds: 

a x {%) =5.2-10- 12 s--H (1) 

with the integration time T in seconds. Equation ((TJ is valid 
for a single satellite pass over a ground station (for integration 
times T lower than 300 s). For longer integrations times 

a x (z) = 2.4- 10~ 14 s-t?. (2) 

The temporal Allan deviation can be related to frequency 
instability as expressed by the modified Allan deviation 
Moda y (x) by 

Moda y {x)= ^ ax ^\ (3) 

We take (Q]) and (f2| as our upper limits for the calculation 
of all perturbing effects in the following sections, together with 
the overall accuracy requirement (maximum allowed frequency 
bias) of 1 • 10~ 16 in relative frequency. 

3. Relativistic model for clocks and time transfer of ACES 

In a general relativistic framework each clock produces its own 
local proper time, in our case X 8 and X s for the ground and 
space clocks respectively. 

In order to model signal propagation between the ground 
and the space stations, we use a non-rotating geocentric space- 
time coordinate system. Thus t = xq/c is the geocentric coor- 
dinate time, x = (xi,X2,x?,) are the spatial coordinates, where 
c is the speed of light in vacuum (c = 299792458 m.s -1 ). We 
denote U(t , x ) as the total Newtonian potential at the coordi- 
nate time t and the po sition x with the convention that U > 
(Sof fel M. et al. 20031) . In these coordinates, the metric is given 
by an approximate solution of Einstein's equations valid for 

I! v 2 

low velocity and potential << 1 and ^ << 1): 

ds 2 = -(l —1 — -cW+ H — -)dx 2 , (4) 

where higher order terms can be neglected for our purposes 
dWolf P. and Petit G. 1995b . 

In this system, each emission or reception event (at the an- 
tenna phase center) is identified by its coordinate time f; (figure 
[T]) and a coordinate time interval is defined by Tjj = tj — ti. We 
define x g , v g and a g respectively as the position, the veloc- 
ity and the acceleration of the ground station, and x s , v s and 
a s respectively as the position, the velocity and the accelera- 
tion of the space station. 
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Fig. 1. MWL principle 



The ACES mission uses two different antennas: one Ku- 
band antenna for uplink and downlink, and one S-band antenna 
for downwards signal only. The antennae are characterized by 
their phase pattern, which describes the position of the antenna 
phase center as a function of the direction of the incoming sig- 
nal. The phase patterns of the MWL antennas have been mea- 
sured in the laboratory at all three frequencies. For example, 
the phase variations with direction of the Ku-band antenna can 
reach up to ^0.1 rad (^1 ps in time). Those variations can be 
corrected from the known phase pattern and a knowledge of 
ISS attitude (ie. direction of the incoming signal). Once cor- 
rected the remaining errors, mostly due to uncertainties in ISS 
attitude, are below the MWL specifications. 

The f\ frequency signal is emitted by the ground station at 
the coordinate time t\ and received by the space station at tj.. 
The f2 and fa frequency signals are emitted from the space sta- 
tion at ?3 and ?5, and received at the ground station at f 4 and fg. 
The third frequency is added to measure the TEC in the iono- 
sphere which allows the correction of the ionospheric delay. 

The MWL is characterized by its continuous way of emis- 
sion. It measures the time offsets between the locally generated 
signal and the received one. It provides three measurements (or 
observables) of the code (one on the space station, two on the 
ground) and three measurements of the phase of the carrier fre- 
quency at a sampling rate of one Hertz. 

An observable is related to the phase comparison between a 
signal derived from the local oscillator and the received signal, 
corrected for the frequency difference mainly due to the first 
order Doppler effect (see Bahder T. B. 2003 for details of a 
similar procedure used in GPS). If we consider a particular bit 
of the signal which is produced locally at x p and received at 
X a , an observable is a measurement of the local proper time 
interval between these two events. The actual measurement of 
the phase difference A<t>(T a , x a ) occurs at the single space- 
time point (T a , x a) and is labeled with the arrival proper time 
X a . It can be expressed as 



At(t , ** fl ) = 



A<t>(x a ,l? a ) 



(O 



St = T p - % a 



(5) 



where co is a conversion factor from phase to proper time 
depending on the nominal carrier and code frequencies, and 
8x represents the measurement errors (difference between the 
clock and ideal proper time, measurement phase noise, etc.). 



Considering the experimental uncertainties of the ACES 
mission (see equation (fTJ and (0), we will neglect any terms 
in the relativistic model that, when maximised, lead to cor- 
rections of less than 3 • 10~ 13 s in time. Numerical applica- 
tions are necessary to evaluate which terms can be neglected. 
For this purpose, we consider the International Space Station 
with a circular trajectory at an altitude of 400 km in a plane 
inclined by 51.6" with respect to the equatorial plane. It has 
a velocity v s = 7.7 ■ 10 3 m.s~ l in a gravitational potential of 
U s /c 2 = 6.5 • 10~ 10 . The ground station has a velocity v g = 465 
m.s~ l at a gravitational potential of U g /c 2 = 6.9 • 10~ . 

The tropospheric delay A tropo is considered as independent 
of the frequency of the signal, with a slow variation with time 
and an amplitude of less than 100 ns. 

We assume the density of electrons in the ionosphere N e is 
less than 2 • 10 12 electrons/m 3 . The ionospheric delay A' ono is 
maximum for the frequency signal with an amplitude of less 
than 10 ns. 

The relation between the proper time x and the coordinate 
time t is given to sufficient accuracy by 



dx ( U{t,t) v 2 (t) 

dt \ c 2 + 2c 2 



+ 0(c- 4 ) 



(6) 



Higher order terms of equation © have negligible effects 
at the projected unc ertainty of 1 • 10~ 16 in rel ative frequency 
of the ACES clocks (I Wolf P. and Petit G. 19951) . Note however, 
that some care has to be taken when evaluating the Newto- 
nian potential U(t, x) in © for the ground or space clock 



(IWolfP. and Petit G. 1995h . 

The ACES mission aims at obtaining the variation of the 
desynchronisation between ground and space clocks with time, 
that is to say, the function X s (t) — x s (t). It is evaluated by com- 
bining the measurements performed on the ground and onboard 
the space station and a precise calculation of the signal propa- 
gation times. In order to be able to control T23 (see below), we 
combine two measurements Ax s (x s (t2)) and Ax g (x s (t4)) but 
with x s (t?) ^ T g fa). Then the expression of desynchronisation 
reads (see (Duc hayne L. 2 008) for a detailed derivation) 



x*(t a )-T\t a )= ] -(A: : ■ \: : . 



h c2 



2c 2 



)dt 



(7) 



where t a = and where Ax s (x s {t 2 )) and AT*(T*(f 4 )) 
are the observables respectively from the ground and onboard 
the satellite at the coordinate times t% and f 4 , and where we have 
neglected non-linearities of X s (t) and X s (t) over the interval 
?4 — ?2 (f ew milliseconds). The integral terms result from proper 
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time to coordinate time transformations. They are small correc- 
tions of order 10~ 12 s to the desynchronisation % g {t a ) — T s (t a ). 

However, it is the derivative with respect to the coordinate 
time t of the relation © which has to be studied for applica- 
tions such as the gravitational redshift test or geodesy. Actu- 
ally it has to be compared with the next relation obtained from 
equation ©: 



dz g <ir 1 . 

— - — = -s-[U(tX)-U(t,%) 



dt dt 



v?(0 #0 



(8) 



+ 0(c" 4 ) 



We note that any constant term appearing in the desynchro- 
nisation expression (|7]i will have no effect on the final result ^ 
because of the derivation. 

In (|2l) the difference T\ 2 — T$ 4 needs to be calculated from 
the knowledge of satellite and ground positions and veloci- 
ties (orbit determination). For example, T\ 2 is the time interval 
elapsed between emission from the ground station and recep- 
tion by the satellite of the f\ frequency signal. It can be written 



#12 , 2GM E ( x g (h)+x s {t 2 )+R n 
x g {h)+x s (t 2 )-Ri2 



■In 



■A 1 ™ p0 - 



a iono 
A 12 



+ 0(-), 

cr 



(9) 



where R\ 2 — \\R\i\\ — \\xl{ti) — *g(fi)||> where the logarith- 
mic term represents the Shapiro time delay (IShapiro 1. 1. 1 964) 
(see e.g. Blanchet L. et al. 2001 for a detailed derivation) and 



where A' 



tropo 
12 



and Af 2 "° 



are respectively the tropospheric and 
ionospheric delays on the signal path. 

Only one of the clocks (here we assume the ground clock) 
has a known relation with coordinate time f, so x s (t 2 ) cannot 
be directly obtained from the orbit determination, only x s (ti) 
is known. Relation (O is then modified to 



T 12 = 



D{t x ) D{ tl )X{h) 

2 

c c 
D(h) ( 



2c 3 \ 



U{ h )\\ 2 + D{h).^{h) 



Dt, 



(10) 



2GM E fx g (t i ) + x s (h)+D(t x ) 



x g (h)+x s (ti)-D(h) 
i: ^0(c- 4 ), 



-A\ 2 no + A' mpo 1 



where D(t) = JtJ (t) - x^(t), D(t) = ||Z)(f)||, and where the 
unknown position ~xl{t 2 ) in (O was expanded in a Taylor se- 
ries around the corresponding known position x s (f i ) . The iono- 
spheric and tropospheric terms are related to the signal paths. 



They need to be calculated for R (t), but only D (t ) is known. 
The corresponding correction is related to the motion of the 
space station during the atmospheric delay of the signal (the 
atmospheric equivalent of the "Sagnac" type term (2 nd term in 
([Toll). The corresponding corrections to the tropospheric term 
and to the ionospheric terms at frequencies f\ and f 2 are negli- 
gible but included for the signal (see equation (fT4l). 

Assuming that T 2 ^ < 10~ 3 s , the expression of T\ 2 — T34 
where all terms greater than 3 • 10~ 13 s have been included 
when maximized for ACES can be evaluated by : 



l\2— •'34 = ^ 9 1 TZ77~\ / 23 



c-D{t 4 ) 



D(t 4 ) ( 



_ ^ . \Av(t 4 ).v* s {t 4 ) 

-D(f 4 ).^(f 4 ) + |jA^(f 4 )|| 

7,3 /— > _> 
+ \Av(u).v s (t 4 ) 

-D(/ 4 ).^(f 4 )+2||A^(f 4 )|| 



2D(f 4 ).Aa(f 4 )- ( 



D(t 4 ).Av(t 4 ) 
D(t 4 ) 



(11) 



T ^- I \\Mt 4 )\\ 2 -D(t 4 ).Aa(t 4 ) 



2cD(t 4 ) 



D(t 4 )Av{t 4 ) 



v D(t 4 ) 

1 a iono a iono 



where Av(t) = ~v 8 (t) — ~v s (t) and Aa(t) = ~a ' g (t) — ~~a s {t). 

To obtain (fTTT i we have applied ( fTOb for the upward and 
downward signals and expanded all positions and velocities 
in Taylor series around their values at f 4 , which can be ob- 
tained from the time of measurement ("label" of the observable 
At* (t 8 (t 4 )) in ©), on the ground. 

The difference T\ 2 — T^ 4 of upward and downward signals 
at /1 and f 2 allows to eliminate to first order delaying and re- 
straining factors such as range (D/c), troposphere or Shapiro 
effects. Due to the asymmetry of the paths, that cancellation 
is not perfect, and there are some terms left (equation (fTTI) ) 
which depend on orbit determination as well as on the coordi- 
nate time interval T 2 j elapsed between reception and emission 
at the phase centre of the MWL antenna onboard the ISS. 

The time interval T 2 ^ can be controlled by "shifting" the 
time of the measurements on the ground or on board the space 
station (e.g. shifting T s {t 2 ) in Ax s (T s (t 2 )) of equation (|7]i) when 
post-processing the data. That allows the control of the differ- 
ence between emission and reception at the clock. But 723 in 
(fTTT i is defined at the antenna phase centre, therefore its con- 
trol requires the knowledge (calibration and in situ measure- 
ment) of the instrumental delays (cables, electronics, etc..) in 
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the MWL space segment. So the control and knowledge (uncer- 
tainty 8T23) of is determined by the uncertainty of the cal- 
ibration and measurement of internal delays, which will play a 
role in the following. 

The ionospheric effect A lo no on a radio signal of frequenc y 
/ can be modelled as follows ( Bassiri S. and Haii G. A. 1993b : 



40.3082 



c-f 2 
7527 



N e dL 



P 



(B . k)-N e dL + 0( 



1 



(12) 



where Bo is the Earth's magnetic field (||fio|| ~ 3.12- 
10~ 5 T) and k is the direction of propagation of the signal. 

The electron density integrated along the signal trajectory, 
/ N e dL, is usually referred as the Total Electron Content (TEC). 
The difference of ionospheric delays appearing in equation ( ITTb 
is then 



M2 



\iono / 



1 ,40.3082 „ 
- T ) TEC 

fl c 



7527 f B .D(t 2 ) 



fl 
7527 



D(t 2 ) 
D{U) 



N e dL 



N,dL- 



(13) 



The second and third term in (Tl3T > represent the effect of the 
Earth's magnetic field and amount to at most 0.5 ps for ACES, 
which is almost negligible. Therefore, an only very rough esti- 
mation of these terms is sufficient. 

The TEC is determined by combining the two downwards 
signal observables. The time interval Z46 elapsed between the 
receptions of the two signals depends on the differences of their 
internal delays at the emission from the space station (i.e. 
and on the difference of their propagation time. Assuming that 

< 1 H s > the expression of the Total Electron Content is then 
given by the following expression : 



TEC 



1 



40.3082 fl-fl 



D(t 4 ) 



fifl 40 - 3 

AT A '(T A '(f 4 ))-AT«(T S '(f6)) 

D{t A )A-y(t A ) T 46 \ 
D{U) 
D(t 4 ).v* s (t 4 ) 



D(t 4 )-c 



c) 
1 



/ I 



(14) 



0{ 



The last term in (fl4l is negligible when used for the time 
transfer (when inserted into ( fT3l ) but may be significant for the 
study of the TEC itself. 



In summary, a reliable orbit determination is required for 
two main reasons. On one hand to calculate precisely the cor- 
rections in equations ( fTTT i. On the other hand, to evaluate cor- 
rectly the terms on the right hand side of equation dHJ ie. the 
second order Doppler and gravitational redshifts. 

In addition, we also need a precise knowledge of the time 
interval T23, (ie. of the onboard internal delays) in order to be 
able to calculate the corresponding terms in ( fTTT i with sufficient 
accuracy. 

The aim of this work is to estimate using simple orbit error 
models, which levels of uncertainty on orbit determination and 
calibration of internal delays (knowledge of T23) are required 
to reach the required performances. For that purpose only the 
leading terms in ( fTTT i are required ie. 



'12' 



734 



D(t 4 ).v^(t 4 ) 



D(t A )Av{H) 



T23 



(15) 



c-D{t 4 ) 

which, together with equation ^ for the relativistic frequency 
shift, is sufficient to derive the maximum allowed uncertain- 
ties on orbit determination and internal delays in order to stay 
below the limits given by (UJ and (0. 

4. Effects of orbit errors on clock comparison 

We now use equations ( fT3T > and (JHJ to express the effects of 
station trajectory and time calibration uncertainties on the time 
transfer and on the relativistic frequency shift. In this section 
we make no assumptions on orbit determination errors, our re- 
sults being completely general and valid up to the neglected 
terms as described below. 

We note (X a (t), ~V a (t)) and (X' a (t), V' a (t)) respectively 
the true and computed trajectories of the antenna phase center, 
and (X c (t), V c {t)) and (X' c .(t), V' c {t)) respectively the true 
and computed trajectories of the clock reference point. We also 
define (X a (t), V (t)) the true trajectory and the true veloc- 
ity of the center of mass of the ISS. These five trajectories are 
expressed in the non-rotating geoce ntric frame (GCRS, Geo- 
centric Celestial Reference System) (ISoffel M. et al. 2 003). 

On one hand, the error in the time transfer is related with 
the uncertainties of T\2 — T^ 4 , and can be obtained from the sim- 
plified equation < fT3T >- It is then dependent on the ground and 
space station trajectory knowledge, on the value of 723 and on 
the uncertainty on this parameter bTn- As said before, a pre- 
cise knowledge of the time interval 723 is related to the internal 
delay calibrations. The error on 71 2 — 734 is 



8D.vt+ D.dvt 
8 (Z12 - 734) = 2 5-= 8 - 



D.Av 
c-D 



8T7 



2j 



8D.AV D.SA~V 



c-D 

D .A~v 8D 
c-D D~ 



c-D 



(16) 



723- 
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In general ground station uncertainties are below 10 cm. 
Thus the uncertainty on ground station position is negligible 
with respect to the ISS position errors and the knowledge of 
the vector D is only related to the uncertainty of the space 
station reference point which is the antenna phase center. Then 

we have 8 D = X a X' a . 

The previous equation can then be written as : 



5(r 12 -r 34 ) = 2^ i 



D.Av 



8T 7 



23 



X a X' a Av 
c-D 



c-D 
D dXgX' a 
c-D' dt 



(17) 



DAv\\X a X' a \ 
c-D D 



?23- 



We note that equation (fTTT i depends the antenna phase cen- 

> 

ter position error {X a X'^, the coordinate time interval T23 and 
the uncertainty on this time interval 5723. 

On the other hand, the clock relativistic correction along 
a trajectory is obtained from equation ©. It depends on the 
position and the velocity of the reference point, in this case 
the clock. We need to express the error on the reference point 
frequency - ie. the frequency difference between the true clock 
position and the computed clock position - in order to compare 
its Modified Allan deviation with the specifications. 

The gravitational potential c an be evaluated on a give n tra- 
jectory with sufficient precision dWolf P. and Petit G. 1995b us- 
ing gravity models (eg. GRIM5 or EGM96). The error on the 
frequency shift at the clock position is given by : 



nidi, ( dz dT 



1 



u{t,x c )-u{tX) 



(18) 



Vi - V'- 



The frequency difference between the reference point posi- 

— > — > 

tion X and the ISS center of mass X is given by : 



,dx . ,dx „ 1 . 

^h-(^h =--2[ u (^ x )- u ^ x o) 



dt 



V z -V, 



(19) 



The trajectory of X is the solution of the differential equa- 



tion 



d z X t 
dt 2 



= r P 



(20) 



where r p is the acceleration due to the gravitational poten- 
tial and r s is the acceleration due to other effects (e.g. surface 
forces like air drag and radiation pressure). 



Using 



and to first order 



r p = Grad(U), 



U(t, X ) = U(t, X ) + T P (X ) XoX, 



(21) 



(22) 



and multiplying equation ( f20l > by the vector X„X (which 
is the position of the clock with respect to the ISS center of 
mass), we can substitute the difference of gravitational poten- 
tial in equation ( fT9l by : 



U(t,X)-U(t,X )= T P {X ).X X 
d 2 X 



(23) 



dt 2 



.XgX— r S-X () x 



Then we obtain : 



(-)v-(-hr 



dV 1 -> dX a X 



1 , dXstX . o t \ 



(24) 



which can be simplified to : 



(-)">-(-)-> 



dt 



+ 2< 



\ f dXoX 2 



dt 



rsXnX 



(25) 



In this expression the term \ {^jjf-) 2 — T s-X„X can be in- 
terpreted as the relativistic correction for the clock referenced 
to the local ISS frame. In fact, for an observer in the ISS frame, 
the non gravitational acceleration F 5 produces an acceleration 
which can be seen by this observer as coming from a "gravita- 
tional potential" — F s-X () X. The term (v o .X X^j is the po- 
sition of the clock with respect to the ISS center of mass pro- 
jected on the ISS velocity direction (along track component). If 
there are no external accelerations, and if the velocity of X in 
the ISS frame remains small, only this term is present. 

Combining the previous equation with the same expressed 

— > 

in X', we obtain : 



8( 



dt >x c 



(y .xjfy+-( 



1 / dX Q X c 2 



dt 



(26) 



In order to simplify equation ( I26l l. we evaluate the order of 
magnitude of the different contributors appearing in this equa- 
tion. 

To investigate the importance of the non gravitational term 
Y s - XcXc > the drag has been modeled along a reference orbit of 
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ISS drag acceleration 




22000 22100 22200 22000 22400 22000 22000 22700 22000 22000 20000 



Fig. 2. Estimation of the non gravitational acceleration of the 
ISS vs time (in days) 

the ISS. A period with important solar activity has been chosen 
in order to evaluate the worst case (see figure[2]i. 

To estimate its effect on formula (1261 , the acceleration has 
been multiplied by a 10 meter bias, by a 10 meter random 
noise or by a 10 meter sinusoidal function at orbital period, 
corresponding to possible attitude and orbit error effects of 
the ISS. The Allan deviation stays below 10 -21 , which is to- 
tally negligible here. Also, this term has no effect on the fre- 
quency accuracy at the 10~ 16 level aimed at by ACES. In 
addition, the residual term of the second order Doppler shift 



( dA g A c \ 
V (it > 



must be computed with the GCRS 



l 

trajectories. The order of magnitude of these terms can be eval- 
uated as ~ (aSaQ. 2 /c 2 ) with £2 the orbital angular frequency, 
a = \ \X X C \\, and 8a its error. For a < 100 m and 8a < 10 m 
this effect is also totally negligible. 

The only important term for the performance evaluation is 

thus the along track term 
can be written : 



V .X C X' C j , and equation 



,dx 



,dx , 



dt ( dt ^' 



(27) 



So only the component of the clock position error X C X' C pro- 
jected on the velocity of the ISS V plays a role. This can be un- 
derstood considering for example a purely positive radial com- 
ponent. In this case we underestimate the gravitational potential 
but overestimate the velocity, and the two cancel. We underline 
again the fact that the derivations in this chapter are valid for 
any type of orbital errors, up to the neglected terms described 
above. 

The scalar products of vectors can be evaluated in a local 
frame : for example it may be useful to study them in the local 

— y — > — y — > 

orbital frame ( R , T , N) defined with R the unit vector be- 
tween the Earth's center and the space station, N orthogonal to 
R and the inertial velocity, and T orthogonal to R and N . 



The ISS attitude is assumed to be roughly constant. The 
clock position error in this frame is represented by a constant 
vector (error in the position of the clock relative to the centre 
of mass), and by perturbations which reflect : 

- attitude uncertainties (rigid body behaviour of ISS). If we 
take ±5 degrees uncertainty and a distance of 10 meters, 
this leads to 0.87 meters amplitude error for the clock posi- 
tion. 

- ISS deformations : they are due to thermoelastic effects 
which will be mainly at orbital period (eclipses and sun 
orientation). We suppose that these effects are below one 
meter amplitude. 

- vibrations : here we suppose that the 'vibrations' are re- 
lated to frequencies higher or equal to the first ISS eigen- 
frequency which is higher than the orbital period. The cor- 
responding displacements remain small and are expected to 
stay below one meter. 

Either of these effects is negligible when inserted into dTTb and 
d27b . So for our purposes the only significant contribution to 
the trajectory errors of the antenna phase center and the clock 
come from the position errors of the ISS center of mass. 

In summary, the errors on the time transfer and on the rel- 
ativistic correction have been expressed as functions of the tra- 
jectory knowledge through equations ( fT71 ) and j27\ . Moreover, 
only the error in the knowledge of the ISS center of mass posi- 
tion has an importance in the relativistic correction. 

5. Orbit error examples 

In this section we describe two simple models for ISS orbit 
errors as examples to investigate numerically the effect they 
have on ACES performance. 

We consider two methods used to evaluate the orbit of a 
satellite. On the one hand, the dynamical method takes into 
account the equations of motion to restitute the trajectory. It 
fits the measured data towards an orbit satisfying the principles 
of celestial mechanics and estimates a number of parameters 
which give the best fit trajectory. On the other hand, the kine- 
matic method uses directly the measurements of the satellite 
position (eg. GPS data) with no a priori assumption of the form 
of the orbit. In both cases the resulting orbit errors are in gen- 
eral not easily described by a simple model, but for our pur- 
poses we use two examples that should approximately reflect 
the orders of magnitude of the main orbit error contributions. 

For dynamical orbit determination the differences between 
true and computed trajectories of the ISS center of mass are 
expected to have specific structures. For example an eccentric- 
ity error gives no long term effects, but periodic errors can be 
important and the radial, along track and velocity errors are 
correlated. For weak eccentricity orbits, the difference between 
two real orbits is given by the Hill model (or the Clohessy- 
Wiltshire model) which is an expansion of uncertainties with 



respect to a reference ci rcular orbit dColombo O. L. 1986 



Colom bo O. L. et al. 20041) . If we suppose there exists a set of 
parameters which perfectly describes the true orbit of the ISS, 
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then this model is a first approximation to the errors of the cen- 
ter of mass between the computed and the true orbits. 

For kinematic orbit determination one expects little a priori 
correlation between the different orbit error components. We 
therefore use a simple independent random noise on each com- 
ponent as a first approximation. 

In the Hill model the position uncertainties along radial, 
tangential and normal axis (in the local ISS frame) are given as 
follows: 



radial axis : 8R(t ) = ■ cos(ilf + (p R ) + c R 



each couple (X, T23) gives a deviation staying under the speci- 
fications, and the prohibited area. 



tangential axis : 8T(t) = —X ■ sin(i2f + (p R ) 
normal axis : 8N(t ) = Y ■ cos(i2f + (p^) 



-Q.CR-t 
2 



(28) 



where X, Y, cr and d R are amplitude coefficients, and where 
Q. is the orbital angular frequency. For our purpose, bias (cIr) 
plays no role and the linear term (cr) depends on arc length. 
Basically the longer the observation duration is, the smaller this 
coefficient becomes. The main feature of this error model is to 
take into account the error correlations in the orbital plane. For 
instance, a positive radial bias leads to a negative error on the 
tangential velocity: the satellite is delayed with respect to the 
reference orbit. 

6. Numerical results 

We now use the previously described error models to calculate 
the corresponding constraints imposed by the ACES stability 
(equations ([TJ and (|2]i) and accuracy requirements, via equa- 
tions (T% and 1271 . We first consider the Hill model, followed 
by results for a simple random noise. 

We consider orbit error models in the local (R,T,N) 
frame, so we need to transform them to GCRS and determine 
the uncertainties in this frame on position and on velocity pa- 
rameters. We consider an ephemeris of the ISS correspond- 
ing to the 20 th of May, 2005 and a ground station based in 
Toulouse, France (43°36'/V, 1°26'£). Actually, this station has 
been chosen as the master ground station of the ACES mis- 
sion. All station parameters and their uncertainties have to be 
expressed in the same frame (ie. GCRS). 

We first consider the error equation ( TP71 ) on the time trans- 

> 

fer. We choose the signs of the independent parameters (X a X' r 
723 an d 8T23) so as to maximize the resulting temporal Allan 
deviation. The calculated deviation has to be compared with 
the MWL specifications. 

Assuming we have no error on T23 (ie. 8T23 =0 s), for all 
values of factor X (or Y) of equation (l28l l. it is possible to de- 
termine the maximum value of the time interval 723 for which 
the temporal Allan deviation remains under the specifications. 
With numerous values of X, we calculate a bound which marks 
out two different areas : the allowed uncertainties area in which 




Factor X (m) 



Fig. 3. Maximum allowed value of T23 as function of the scale 
factor X to comply with the specifications, assuming 5723 =0 

Figure [3] shows that, the smaller the time interval 723, the 
larger the allowed uncertainty on the space station position. 
This result provides a way to combine upwards and downwards 
signals to allow the maximum uncertainty on space station po- 
sition in order to comply with the specifications. The most fa- 
vorable situation is when the reception at the antenna phase 
center of the space station corresponds to the emission at the 
same place ie. t2 = ?3 (see figure HJ. This way of combining 
signals is named the "A configuration". To work with param- 
eters in the asymptotic area (see figure O requires 723 to be 
below 10~ 6 s. 



Space station path 



k=fc 




Ground station 
^trajectory 



EARTH 

Fig. 4. The "A configuration" is the way to combine upwards 
and downwards signals which allows the maximum uncertain- 
ties on the space station orbit determination 



Then if we plot the maximum value of 8T23 for all values 
of the factor X, there will appear two asymptotic values we 
cannot cross if we want to stay within the specifications (see 
figure [5]). Basically a compromise between the knowledge of 
the space station trajectory and the precision of the internal de- 
lays calibration must be achieved owing to the maximization 
of the Allan deviation. We will evaluate the maximum allowed 
errors on these two parameters if no other errors are present. 

We search for the asymptotic value of factors X and Y 
which comply with the specifications for all phases (<p^, <pw) 
when we have no error on 723 . The asymptotic value for orbit 
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Fig. 5. Maximum allowed value of 8T23 as function of the scale 
factor X to comply with the specifications, assuming T23 = 
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Fig. 7. Modified Allan deviations of the redshift error forX=14, 
16, 18, 20 m 



determination is obtained forX, Y = 4200 m which corresponds 
to a 4.2 km error on the tangential and normal axes, and to a 
2. 1 km error on the radial axis. 

The asymptotic value of the time calibration does not de- 
pend on orbit determination uncertainties. So it is independent 
of the phases <Pr and (p^- Then we can draw the temporal Allan 
deviations for different values of 5723 (see figure [6]). We find 
that 8T23 must stay below 1 .06 • 10~ 7 s. 



Time Stability 




Time (s) 



Fig. 6. Temporal Allan deviations for X = 0, T23 = s and 
ST 23 = [102, 104, 106, 108] ns 



The requirements for several passes have also been inves- 
tigated. In this case, the calculated deviations have to be com- 
pared to the specifications given by (0. The results showed 
that the requirements on orbit determination and time calibra- 
tion are less stringent for several passes than for a single pass. 
Therefore if specifications are respected for a single pass, spec- 
ifications for longer integration times are also respected. 

Now we evaluate requirements on orbit determination con- 
sidering the relativistic frequency shift. We search for the max- 
imum value of X to comply with the specifications (fl3 and (|2j. 
Equation (|27T i is evaluated with the error model d28l . and its 
Allan deviation is calculated for different values of X. For in- 
tegration times greater than one thousand seconds, these Allan 
deviations are independent of the phases (pg and q>pj. 



Figure [7] shows that, if the factor X in ( |28l is equal to 16 
m ie. if we have an eight meter error on the radial axis and 
sixteen meter error on the tangential axis, then we comply with 
the specifications. 

Because of the projection of the position error along the ISS 
center of mass velocity (see equation (l27l>). the requirement on 
the factor Y is orders of magnitude less stringent. The bound 
is given by the asymptotic value coming from the time transfer 
(see equation ( TPTT i and figure|5]l. 

We now turn to the second example of orbit determination 
errors: independent random errors on R, T and N. We con- 
sider white noise of amplitude (standard deviation) 8R, 8T, 8N 
at one second sampling intervals on each of the three compo- 
nents. 

Similarly to the Hill model, the by far more stringent con- 
straints come from the effect of these errors on the relativistic 
frequency shift via equation (l27l i. In particular, a white noise 
8T in equation (l27T i will translate into a temporal Allan vari- 
ance given by 



VqST 



(29) 



which satisfies the requirements ([T|i and (O for all integration 
times when 8T < 60 m. 

Because of the projection of the error on Vq in d27b the rel- 
ativistic frequency shift imposes virtually no limits on the nor- 
mal and radial components of the errors. Those are limited by 
the effect on the time transfer given by equation ( fTTl ). Assum- 
ing < 10~ 6 s and 5723 < 10~ 7 s (see above) the first term of 
that equation is by far dominant. Maximizing the dot product 
in that term we obtain an upper limit of 1 .4 km on 8R and 8N 
in order to stay within specifications for all integration times. 

Finally, we consider the accuracy requirement of ACES i.e. 
10~ 16 in relative frequency when averaged over ten days. From 
the integral of equation ( f2Tb this implies that the tangential 
component of the position error (X C X' C in d27i i) cumulated over 
ten days needs to remain below one kilometer (including for 
example the linear term along the tangential axis in d28l)). This 
is unlikely to raise any difficulty, if the much more stringent re- 
quirements from periodic or random errors (see above) are met. 
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We note that the accuracy requirement is related to a require- 
ment on the knowledge and conservation of the total energy of 
the orbit. An error in the total energy will show up as a velocity 
bias, and thus as a linear term along the tangential axis. The 
non-conservation of the total energy (non-gravitational accel- 
erations) was found to be negligible in section @] Nonetheless 
possible long-term effects, if they exist, will lead to a slow vari- 
ation of the tangential error. The constraint on both, the error in 
the total energy and the effects due to energy non-conservation 
is provided by the accuracy requirement, i.e. that the tangen- 
tial error cumulated over ten days needs to remain below one 
kilometer. 

7. Conclusion 

We have derived detailed and general (within the specified as- 
sumptions) expressions for the calculation of the effect of or- 
bit determination and instrumental calibration errors on the 
stability and accuracy of ground to space clock comparisons. 
These expressions were then used on the specific example of 
the ACES mission to derive maximum allowed errors given the 
stability requirements of that mission. For that purpose we have 
used two simplified orbit error models (Hill model, and random 
noise), setting limits on the appropriate parameters of those 
models. Although neither of those models is likely to correctly 
reflect all of the ISS orbit errors, they are nonetheless expected 
to provide correct orders of magnitude for the maximum al- 
lowed orbit determination uncertainties. They are summarized 
in table Q] where in each column we have provided the more 
stringent result from the two models. 

Furthermore, we have identified an optimal way to combine 
upwards and downwards signals (the A configuration), which 
allows for the maximum orbit determination errors. This then 
provides a constraint on the maximum allowed value of T23 < 
10~ 6 s and its maximum allowed uncertainty of 8T23 < 10 -7 s, 
which constrains the required knowledge of the sum (transmis- 
sion + reception) of instrumental delays between the antenna 
phase centre and the onboard clock (also provided in table[T]i. 

Table 1. Requirements on orbit determination and internal de- 
lay calibration 



SR/m 


ST /m 


SN/m 


5r 23 /s 


8 


16 


1400 
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Finally, the accuracy requirement of ACES was used to 
constrain long term (1=3 10 days) linear drifts and other slow 
variations of the orbit determination errors. We find that the to- 
tal tangential error cumulated over 10 days needs to stay below 
1 km to comply with the accuracy specification (10~ 16 in rel- 
ative frequency after ten days averaging), which should pose 
no particular problems given the more stringent short term re- 
quirements (see table[T|i. 

In conclusion the requirements on orbit determination are 
significantly less stringent than the initial 'naive' estimate (one 



meter error for 10 in relative frequency), which is mainly 
due to partial cancelation between the gravitational redshift and 
the second order Doppler effect in the relativistic frequency 
correction of the onboard clock. 
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